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Abstract 

A new phenomenon, entrainment of chaos, which is understood as a seizure of an irregular behavior by limit 
cycles, is discussed. As a result, chaotic cycles appear if the chaos amplitude is small. Otherwise, the chaos is 
| not necessarily cyclic, if the perturbations are strong and/or diameter of the limit cycle is small. Sensitivity 
as a main and a unique ingredient is considered and, in addition, period-doubling route to chaos is chosen for 
extension. The results may be of strong importance for engineering sciences, brainwaves and biomusicology 
phenomena as well as can be developed for hydrodynamics. Theoretical results are supported by simulations 
and discussions over Chua's oscillators, entrainment of chaos by toroidal attractors and controlling problems. 
Moreover, through an example, by means of the Lyapunov functions method, a chaotic attractor is provided. 

Keywords: Orbital stability, Limit cycle, Chaotic attractor, Sensitivity, Period-doubling cascade, Hopf 
bifurcation, Toroidal attractor, Chua's oscillator, Chaos control 



1 Introduction 

It was Christiaan Huygens who observed that two of his pendulum clocks mounted next to each other on the 
same support often become synchronized. He called this synchronization tendency as "entrainment". One 
can say also about the practice of entraining one's brainwaves to a desired frequency, that is, brainwaves 
entrainment [TJ [5] or biomusicology as the synchronization of organisms to an external rhythm. Entrainment 
phenomenon is known also in hydrodynamics as the movement of one fluid by another. Results of our paper 
ensure us to say about entrainment of chaos, which can be applied, for example, in economics theory to 
achieve irregularity in business cycles [5] , and to obtain chaotic cycles in electrical circuits by means of Van 
der Pol equations j4j and in chemical oscillators such as Belousov-Zhabotinsky reactions in continuous-flow, 
stirred tank reactors [5]. Since, starting with experiments of Huygens, the phenomena is mainly related to 
cyclic motions, more precisely, entrainment of chaos by limit cycles is under discussion in our paper. Thus, we 
demonstrate that the idea of entrainment is not concerning only frequency, period or phase OH], but chaos, 
also. In our study the concept of entrainment is less symmetric than that for Huygens synchronization. 
It is, rather similar to the phenomenon in biomusicology or hydrodynamics. Nevertheless, unidirectional 
couplings have been intensively considered in physical studies [7]-|llj. We hope that our investigations will 
be useful first of all for neural sciences, then they are of some interest for chemistry, mechanics, electronics 
and population dynamics. Shortly speaking, they can be applied and developed for any research, where limit 
cycles have been observed. 

In general, we consider systems, which admit well arranged, steady motions and subject to them external 
irregular perturbations. The question is, how their motions are tuned to the external influences? If the 
external forces are regular, then reactions are well analyzed in classical literature. For periodic perturbations 
- periodic responses, if perturbations are almost periodic, then responses are almost periodic, etc |12| 113]. 
What about chaotic perturbations? An answer to this question has been introduced in our papers |14]-|20j. 
where we discuss the problem if systems are with asymptotically stable equilibriums, and found that solutions 
admit the same type of chaos as perturbations do. 
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The paper [14) . where we discuss the mechanism of morphogenesis of chaos, informs us about replication 
of specific types of chaos, such as Devaney and Li-Yorke chaos and chaos obtained through period-doubling 
cascade, in continuous-time systems with arbitrary high dimensions. In this process, we take into account 
the generator-replicator systems such that the generator is considered as a system of the form 

x' = F(t,x), (1.1) 

where F : R x R m — > R m is a continuous function in all its arguments and the replicator is assumed to have 
the form 

y' =Ay + g(x,y), (1.2) 

where g : R m x R™ — > R™ is a continuous function in all its arguments and the constant n x n real valued 
matrix A has real parts of eigenvalues all negative. 

The rigorous results of the morphogenesis (extension) mechanism emphasize that system (|1.2[) is chaotic 
in the same way as system (|l.ljl . Replication of chaos through intermittency is also demonstrated through 
simulations in paper |14| . where one can find new definitions for chaotic sets of functions, and precise de- 
scriptions for the ingredients of Devaney and Li-Yorke chaos in continuous-time dynamics, which are used as 
tools for the study of morphogenesis of chaos. 

Distinctively from the morphogenesis mechanism where we considered chaos near fixed points, in the 
present article, we take into account systems with orbitally stable limit cycles, and perturb them with chaos. 
As a result we obtain chaotic cycles, that is, motions which behave cyclically and chaotically simultaneously. 
Since one achieves appearance of chaotic behavior near limit cycles, the sizes of the attained chaotic attractors 
are possibly larger than the ones considered in paper [14] . Even one can say that the chaos extension procedure 
developed in the present article is global, while it is local in the former one. 

The theoretical novelty of our analysis lies in the fact that we apply the complex technique of the proof 
of Andronov-Witt theorem [3T] to indicate the boundedness of solutions, and describe in detail the cyclical 
behavior of the trajectories of the chaotically forced system in a neighborhood of the limit cycle. Extension 
of sensitivity and chaos through period-doubling cascade are also elaborated rigorously in the paper. 

Entrainment of chaos should be understood as the seize of chaos by the vector field near the limit cycle. 
A novelty of the entrainment process, which is important for applications, is the combination of a chaotic 
behavior with a cyclical one. Cyclical behavior of the motions can disappear if the applied perturbations are 
strong and/or diameter of the limit cycle is small. This type of behavior is discussed in the third section of 
the paper. We prove presence of chaos through indicating the period-doubling cascade, and sensitivity, and 
recognize that observation of other ingredients of chaos, transitivity, density of periodic motions occurs as 
a more difficult task, than we presupposed. Nevertheless, it is known that [22j [23] sensitivity is the main 
ingredient of chaos, since it assumes unpredictability. We leave the discussion for the attractiveness of chaotic 
solutions, but we consider the problem in an example. 

It is very important to comprehend how one can develop present results in further investigations. From this 
point of view, in the sixth section of our paper, we place simulations of chaotically perturbed toroidal attractor, 
which results as a chaotic torus, but we do not give a rigorous mathematical basis for this phenomena. 
Moreover, we perform our method to perturb chaotically stable Chua's oscillators to attain a new chaotic 
Chua's attractor. Finally, to observe unstable periodic solutions present in the generated chaos, we give an 
example for the application of Pyragas control method. 

To present the main idea of the paper, we continue with an example of an oscillating chemical reaction as a 
possible model for our theory. Paper [M] considers the chlorine dioxidc-iodinc-malonic acid (CIO2 — 12 — MA) 
chemical reaction which possesses the following three component reactions: 

MA + I 2 — > IMA + I- + H+, 

CIO2 + 1' — > C70 2 - + \I 2 , (1.3) 
CIO2 + 4J~ + AH+ — > CI- + 2I 2 + 2H 2 0. 
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The rate equations corresponding to reactions (|1.3p are given by 



d[I 2 ] k Xa [MA)[h] 



dt k lb + [I 2 ] ' 

= -k 2 [Cl0 2 }[I-}, (1.4) 



d[ClQ 2 



dt 

-^—^-k 3a [Cl0 2 ][I }[H } , 

where fc la , kn,, k 2 , k 3a , k 3b are the rate constants and s 1 / 2 represents the level of above which the inhibitory 
effect of iodine ion becomes significant. 

Due to the complicatedness for analytic calculations, after approximating the concentrations of slow 
reactants MA, I 2 and C10 2 as constants and making reasonable simplifications and nondimcnsionalizations, 
Lengycl, Rdbai and Epstein [24] reduced system (|1.4[) to the 2— dimensional system 



u[ = a — ui — 



bu, 1 



4uiu 2 



?<2 



,2 ' 

(1.5) 



l + i 

where u\ and u 2 represent the dimensionless concentrations of I~ and C10 2 ions, and the parameters a > 
and b > depend on the empirical rate constants and on the concentrations assumed for the slow reactants. 

It can be verified that for a given value of the parameter a, system (|1.5p undergoes a Hopf bifurcation at 
the parameter value b = bo = 3a/5 — 25/a such that when b > bo, all trajectories spiral into the stable fixed 
point (u\,u 2 ) = (a/5, 1 + a 2 /25), while for b < bo trajectories are attracted to an orbitally stable limit cycle. 
In that case, if wc consider system (|1.5[) with the coefficient a = 11, Hopf bifurcation occurs for bo = 238/55 
and consequently, for b = 2.1, an orbitally stable limit cycle takes place [25] . 

Next, we take into account the Birkhoff-Shaw chaotic attractor [21)1 |2"T] which is generated by the system 
of differential equations 

x[ = 0.7x2 + 10xi(0.1 - xfj, , . 

x' 2 = -xx + 0.25sin(1.57t). [ ' 

We note that system (|1.6[) is in analogy with Van der Pol type equations such that if the periodic forcing 
term 0.25 sin(1.57i) is removed from the second line of (|1.6|) . then under the transformation x — \/T0x 2 , x' — 
— \/T0xi, one attains the Van der Pol equation 

x" + a (x 2 - \)x' + WqX = Ao sin(fc <), (1-7) 

with the coefficients ao = l,u>o = \/0.7 and Ao = 0. The periodic forcing is used for the velocity term in 
system (|1.6[) , while the acceleration is periodically driven in Van der Pol oscillators in the form of equation 
(|1.7[) . This type of forcing is unusual in mechanical systems, but can be observed in electrical or chemical 
systems [26] . 

Making use of system (|1.6j) as the source of chaos and combining with system (|1.5|) in a unidirectional 
way, we set up the following 4— dimensional system 

x'i =0.7x2 + 10xi (0.1 -xl), 
x' 2 = -xi +0.25sin(1.57i), 

^3 = 11-*,- ^ + 1-2 tan (|), (1.8) 

x' 4 = 2.1x3 ( 1 —^r) +0.8x2. 

V 1 + H J 

Practically, it may not be possible to obtain a chemical reaction which admits the rate equations such 
as in system (|1.8[) . However, in this exemplification, our purpose is to realize the feasibility of extension of 
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irregular behavior in system (|1.5[) from an arbitrary source of chaos or chaotic functions, in the case of a 
unidirectional perturbation is actualized as in system (|1.8[) . 

The rigorous results of the present paper indicate that system (jl.8p possesses chaotic motions in the 
4— dimensional phase space and the projection of the chaotic attractor of system (jl.8p on the X3 — 24 plane 
takes place around the orbitally stable limit cycle of system (|1.5[) with the specified coefficients. The appear- 
ance of chaotic behavior around the limit cycle is an indicator of entrainmcnt of chaos. 

To provide an illustration for the mentioned behavior of system (|1.8[) . in Figure [TJ the 2— dimensional pro- 
jections of the trajectory of system (jl.8p with the initial data xi(0) = 0.2, #2(0) = 0.3,23(0) = 0.75, £4(0) = 
4.82 are pictured. Figure [JJ (a) shows the projection on the x\ — xi plane and this picture represents, in fact, 
the Birkhoff-Shaw chaotic attractor produced by system (|1.6I) . On the other hand, Figure [TJ (b) reveals the 
process of entrainmcnt of chaos such that the chaotic attractor generated around the orbitally stable limit 
cycle of system (|1.5[) . where a = 11 and b = 2.1, is illustrated in the picture. Moreover, the irregular behavior 
in the 23 and 24 coordinates of system (|1.8[) generated in process of time is pictured in Figure [5] 



a b 




X 1 X 3 

Figure 1: The 2— dimensional projections on the Xi — X2 and 23 — 24 planes of the chaotic attractor of system 
(|1.8|) are presented in pictures (a) and (b), respectively. The picture in (a) is, in fact, the Birkhoff-Shaw 
chaotic attractor produced by system (|1.6| . while Figure [TJ (b) represents the chaotic attractor generated 
around the orbitally stable limit cycle of system (|1.5j) . where a = 11 and b = 2.1. The irregular structure 
observed around the limit cycle is a manifestation of entrainment of chaos. 
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Figure 2: The graphs of the 23 and 24 coordinates of system (|1.8|) arc presented. The pictures in (a) and (b) 
both present the irregular behavior of the coordinates. 
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2 Preliminaries 



Throughout the paper R and R+ will denote the set of real numbers and the interval [0, oo), respectively. 
Let us consider the system 

x' = F(t,x), (2.9) 

where F : R + x R m — > R m is a continuous function in all its arguments, and take into account the system 

u' = f(u) (2.10) 

where / : R™ — >• R™ is a continuously diffcrcntiablc function. 

To adapt system (|2.10[) in the entrainment process, we perturb this system through the solutions of (|2.9j) . 
and finally obtain the system of the form 

y' = f(y) + w(x), (2.ii) 

where fi is a nonzero real number and g : R m — > R" is a continuous function. 

We mainly assume that system (|2 .9[) admits a chaotic attractor which is understood as a collection of 
functions such that each element of the set is called a chaotic function. More detailed description of chaotic 
functions is provided in P3]. In this case, there exists a positive real number M such that sup tgR+ ||x(f)|| < M, 
for each chaotic solution x(t) of system (|2.9p . The chaotic functions may be irregular as well as regular 
(periodic) such that skeleton of the attractor consists of infinitely many unstable periodic solutions. This 
phenomena is true for both Devaney's concept of chaos [28l [29] as well as for the chaos obtained through 
period-doubling cascade. 

Our another assumption is that the non-linear autonomous system (|2.10[) possesses a non-constant 
T— periodic solution p(t), for some positive real number T, and we consider system (12.111) in a neighbor- 
hood of the orbit 



7 = {u G R n : u = p(t), te[0,T]}. (2.12) 

It is clear that p'(t) is a non-trivial T— periodic solution of the variational system 

v'=A(t)v, (2.13) 

where A(t) = f'(p(t)) is an n x n real, continuous, T— periodic matrix function, and consequently the number 
1 is a characteristic multiplier of system (|2.13p . 

In what follows, we assume that the number 1 is a simple characteristic multiplier of the variational 
system (|2.13[) and the remaining n—1 characteristic multipliers are in modulus less than one. In this case, by 
Andronov-Witt Theorem [3T], the periodic solution p(t) of system (|2.10[) is asymptotically orbitally stable, 
having the asymptotic phase property. 

The following conditions will be used throughout the paper. 

(Al) There exists a positive real number Lf such that — f(V2)\\ < Lf \\yi — 2/2II , f° r all 2/1,2/2 S R n ; 

(A2) There exists a positive real number L g such that Hg^Xi) — (7(2:2) || > L g \\xi — a?2 1| , for all x\, x^ € R m ; 
(A3) There exist positive real numbers Mp,M g such that sup \\F(t, x)\\ < Mf and sup ||.g(a;)|| < M g . 

We emphasize that for an arbitrary solution u(t) of system (|2.10l) . the function z(t) = u(t) — p(t) is a 
solution of the system 

z' = A(t)z + cp(t,z), (2.14) 

where ip(t,z) = f(p(t) + z)-f(p(t))-f(p(t))z. It is clear that ip(t + T,z) = ip(t,z) and ip(t,0) = <p z (t,0) = 0, 
for all t <G R+, and continuity of /' implies that <p z (t, z) = o(l) as z — > uniformly in t £ R+. 
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Similarly, for any given solution x(t) of system (j2.9p , if <f>(t) is an arbitrary solution of the equation 
l/ = f(y) + / i 5( a; (^)); then w(t) = <f>(t) — p(t) is a solution of the following system 

w' = A(t)w + tp(t, to) + ng(x(t)). (2.15) 

Since we take advantage of the proof of Andronov-Witt Theorem, we continue with a brief explanation 
of the technique used in [3T] . 

According to our assumption that system (|2.13p admits the number 1 as a simple characteristic multiplier 
and the remaining n — 1 characteristic multipliers are in modulus less than one, system (|2.13[) has a real 
fundamental matrix $(i) of the form 

*(*) = *(*) (J jL), (2-16) 

where P(t) is a regular, continuously diffcrentiable T— periodic matrix and B\ is an (n — 1) X (n — 1) matrix 
whose eigenvalues are all with negative real parts. 

In the rest of the paper, without loss of generality, we suppose that the coordinates have been trans- 
lated and rotated such that the periodic solution p(t) of system (|2.10[) satisfies that p(0) = and p'(0) = 
(pi, 0, 0, . . . , 0) for some positive real number pi. 

By the help of the fundamental matrix introduced by equation (|2.16p . we define the function 



G(t, s) = < 



/ \ 
$(t) * _1 00, t>s 

^ ? ^ ' v (2 17) 

1 \ 10 
-m[ t<s, 
\ 0„_i / 



where I n -i and O n -\ denote the (n — 1) x (n — 1) identity and the zero matrices, respectively. 

Now, let us choose a real number a > such that Re(ctk) < —a, for all k = 1,2, ...,n— 1, where 
ai, «2, • ■ • , an-i denote the eigenvalues of the matrix In this case, there exists a real number K > such 
that ||G(M)II < A"e- Q (*- S ' for t > s and ||G(M)|| < K for t < s. 

Let L v = Since (p z {t,z) = o(l) as z — >• uniformly in t £ K + , there exists a number S(L V ) > 
such that if ||zi|| < S(L V ), \\z2W < S(L V ), then the inequality \\tp(t, Z\) — tp(t, Z2)\\ < L v \\z\ — Z2W is satisfied 
uniformly in t € R+. 

Suppose that a = (0, 02, 03, . . . , a n ) is an n— dimensional vector which is orthogonal to p'(0). It follows 
from equation (|2.16p that there exists a positive real number K\ such that ||$(t)a|| < K\ ||a|| e~ at , for all 

t € R+. Using the method presented in [3T], one can verify that if ||a|| < then a solution z(t, a) of 

2K\ 

system (|2.14p exists on R + and satisfies the relation 

z(t,a) = $(t)a+ G(t,s)(p(s,z(s,a))ds. (2.18) 



Moreover, the solution z(t,a) satisfies the inequality 

\\z(t,a)\\ < 2Ki \\a\\ er at l 2 , (2.19) 

which is valid for t G K + . 

In the sequel, let us denote by C(^Co) the solution of system (|2.10[) which corresponds to the initial 
condition C(O)Co) = Co- 

For any given initial data Co, the solution C(^Co) provides the relation <C(t, Co) = z(t,a) + p(t), where 
z(t,a) is a solution of system f|2 . 14[) with the initial condition z(0, a) = Co- 
According to equation (|2.18p . we have 

poo 

Ca = ^(0)a+ G(p,s)tp(s,z(s,a))ds = P(0)a-h(a), (2.20) 
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where 

h{a) = P{0)( 1 ° ) f P _1 (s)p(s,z(s,a))cfs, (2.21) 

\ O n -l J JO 

aV t \ 

provided that the inequality ||a|| < y is valid. Since the first column of the matrix -P(O) is [1, 0, . . . , 0], 
it can be verified that the function h(a) can be written in the form 

h(a) = (hi(a2,...,a„),0,...,0^ , 

for some continuously diffcrentiablc function hi. It is easy to show that h(a) = o(\\a\\) as a tends to zero. 
Suppose that £o = (Ci , Cf) • • • j O an d Pij are the coordinates of the matrix P(0), where i, j — 1, 2, . . . , n. In 
this case equation (|2.20[) is equivalent to the following system of n equations 

n 

Cl = ^Plj a j - /ll(a2, • ■ • .On), 

J n 2 (2.22) 

i=2 

The last (n — 1) equations presented in (|2.22p form a linear, regular, one-to-one mapping between the last 
n — 1 coordinates of the vectors a and £o- Expressing <22, • • • , a n in terms of £2 > • • • > C« an d substituting these 
values into the first equation presented in (|2.22[) , one can obtain an equation of the form 

n 

Ci° + E *C° - h(r$, C 3 °, • • • , C°) = 0, (2.23) 

i=2 

where q.;, i = 2, . . . , n, are constants and /i is a continuously diffcrentiablc function such that 

mc 2 °,...,o = o^E(c ) 2 ) 7 

Suppose that S C M n is an (n — 1) dimensional, C 1 manifold determined by the equation 

n 

X! + '^2q i x l - h(x 2 ,x 3 , . . .,x n ) = 0. (2.24) 

i=2 

We note that the manifold S is a hypersurface in a neighborhood of the origin, crossing the orbit 7, which is 
defined by equation (|2.12[) . transversally such that for any solution ((t, Co) of system (|2.10[) with (0 G S, we 
have ||C(ijCo) — p(t)\\ tends to zero exponentially as t tends to infinity. 

Let e = — r — _ . Since the function h(a) defined by equation (|2.2ip satisfies the property h(a) = 
(0)11 

o(||a|j) as a tends to zero, it is possible in this case to find a real number 5(e) > such that if ||aj| < 
min ^5(L V )/(2K\), 5 (e) j- , then the inequality 

h(a) <e\\a\\ (2.25) 

is valid. In the remaining parts of the paper, we suppose that ||a|| < min < — — 5 (e) > . 

I 2Kl J 

The next section is devoted for the existence and behavior of bounded solutions of system (|2.1ip . 
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3 Entrainment and boundedness of solutions 



In Lemma (|3.1[) . we detail the behavior of trajectories near the limit cycle, whose orbital stability is ensured 
by Andronov-Witt Theorem, and then show that the perturbed system admits a set of bounded solutions 
near the limit cycle through Theorem 13. II Thus, we shall prepare a carrier for chaotic solutions. 

Lemma 3.1 For each real number I £ (0,1), there exists a natural number uq = i%o(l) such that for any 
solution C(^: Co ) of system (|2.10[) with Co £ S 1 , the inequality ||£(noT, Co)|| < I ||Co|| is satisfied. 

Proof. Let us fix a solution Co) of system (|2.10[) such that the initial vector Co belongs to the surface S. 
Making use of the equation (|2.20[) and the inequality (|2.25[) one can attain that 

llColl > ||P(0)o|| - \\h(a) 

INI 

Hall 



2iip- i (o)ir 

The last inequality implies that ||a|| < 2 ||P _1 (0)|| ||Co|| ■ In this case, by means of (|2.19j) . it can be verified 
that 

\\C(t, Co) -P(t)\\ = \\z(t,a)\\ < 2K\ \\a\\ < AK X \\P- l {Q)\\ ||Co|| e~^\ (3.26) 

for all t £ K+. 

Now, let us fix an arbitrary number I £ (0, 1). It is possible to find a natural number uq = no(l) such that 
the inequality kK x ||P _1 (0)|| e~ aTno/2 < I holds. Making use of (pT2"6"|) we obtain that 

||C(n T,Co) -K«oT)|| < *K X \\P-\Q)\\ ||Co|| e~ aTn °/ 2 < 2 H Co || • 

Since p{n T) = 0, we have \\((n T, Co)|| < I j|Co|| ■ 
The lemma is proved. □ 

In the rest of the paper, for a given solution x(t) of system (|2.9[) , the function r/ x ^(t, 770) will stand for 
the solution of system y' = f{y) + pg(x(t)) satisfying the initial condition 77^(4) (0, 770) = 770. Furthermore, B r 
will denote the open ball centered at the origin with radius r > 0. 

The next theorem signifies not only a boundedness criteria for solutions of system (|2.11j) , but also their 
cyclical behavior. 

Theorem 3.1 There exist a nonzero real number fj,, a positive real number S and a natural number hq such 
that any solution r] x ( t )(t, rjo) of system (|2.11j) . where x(t) is a chaotic solution of system (|2.9[) and 770 € B$, 
is bounded on R+ such that for all t £ R + the inequality 



rh W {t,m)-p(t)\\ < ^ ( e M« n+ 2)T _ ^ + # (M; ( 3. 27) 



holds, where p= max \\p(t)\\ , and H(S, p) = max(<5 + p) \e 2L ' T , 4K X ||P _1 (0)|| \p + (5 + p) e 2LfT ] } . 
te[o,T] " " 

Proof. The essence of the proof is to determine the numbers p 7^ and 6 > such that for any given 
solution r?a;(t)(*!%)i where x(t) is a chaotic solution of system (|2.9p and 770 € B$, there exists a sequence 
{9i} , &i — > 00 as i — > 00, such that ^(t) (i, ?7o) enters that ball infinitely many times at the moments t = 6*,, 
and consequently the collection of all such functions are uniformly bounded. For our purpose, we will make 
use of the solutions of system (|2.10p which have the same initial data with the chosen function rj x ^(t,rj ) 
at the moments 0j. We note that the sequence {8i} , which will be constructed in the proof, depends on the 
function r] x ^(t, 770), that is, depends on both x(t) and the initial data 770 G B$- 

Since the orbit 7 of the periodic solution p(t) of system (|2.10[) intersects the surface S transversally, there 
exists a real number t\ > such that if \\((t, Co) — p(t)\\ < ei for each t £ [0, 2T], then C(t, Co) intersects S at 



some moment t\ £ [0, 2T]. Suppose that a positive number 6 = d(eij is chosen such that S < eie~ 2LfT and 
let an arbitrary Co £ B$ be given. 

The solution ((t, Co) and the periodic solution p(t) = C(£, 0) of system (|2. 10[) satisfy the relations 

C(i,Co) = Co + [ f(C(s,Co))ds 

and 



o 



?(*) = [ f(p(s))ds, 
Jo 

respectively. Therefore we have 

IIC(*,Co)-p(t)|| < IICo||+ f L f \\asXo)-p(s)\\ds. 

Jo 

Implementing Gronwall-Bellman Lemma [SU] to the last inequality and using ||Co| < S, one can see for 
t £ [0, 2T] that the inequality 

||C(t,Co)-p(*)|| <Se 2L f T 

is valid. Therefore, ||C(^Co) — p(*)|| < e i on the time interval [0,2T], and consequently C(i, Co) intersects the 
surface S at some moment ii(Co) S [0, 2T\. In other words, the point Ci = C(^i(Co)> Co) belongs to the surface 
S. It is clear that ||Cij| < R, where R — t\ + p. 

Now, let us fix an arbitrary number I £ (0, 1). By Lemma l3.ll there exists a natural number no = Hq (^) , 
which is independent of C01 such that 

IICKT + MCo), Co)|| = K(noT, Ci)|| < | llCill < «• (3-28) 
Let us take e = [ | S and choose a nonzero number a such that Ijul < r T , — - rr= 

V 2 y P IPI Mg ^L / („ +2)T _ !j 

Fix an arbitrary solution rj x ^(t, rjo), where x(t) is a chaotic solution of system (|2.9p and 770 £ B$ with the 
number 5 as specified above. Herewith, there exists a number ti (770) £ [0, 2T] such that C(*i( 7 7o) ) f?o) belongs 
to the surface S 1 . 

Making use of the integral equations 



Vx(t)(t,Vo) =Vo+ I f(Vx(t)(s,m))ds + \i J g(x(s))ds 

and 



C(t,Vo) =Vo+ f(C(s,r) ))ds 
Jo 

one can obtain for t £ [0, (no + 2)T] that 

\\v X (t)(t,vo) - C(t,m)\\ < \v\M g t+ / \\vx{t)(s,m) - C( s ,m)\\ ds. 

Jo 

Application of Lemma 2.2 |31| to the last inequality implies that 

\\Vx(t)(t,Vo)-C(t,Vo)\\<\^\M g t+\fi\M g L f f se^-^ds 

Jo 

(e L f l - 1) 



L 



J 



< IpI M 9 ( e L f (n +2)T _ 1 

L f 

< e, 
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for all te [0,(n + 2)T]. 

In this case, we achieve the inequality 

\\Vx(t)(noT + h(r]o),rio) - ((n T + £1(770), Vo)\\ < e , 

and hence by means of (|3.28[) we get 

\\Vx(t)(noT + h(ri ),rio)\\ < \\-q x ( t ){n^T + £1(770), 770) - C(n T + £1(770), rjo)\\ + \\((n T + ii(?7o), Vo)\\ 
<e + lS 
< S. 

In other words, 771 = r] x M (0\, 770) G B$, where Q\ = uqT + £1(770). We note that the point 771 depends on both 
the initial data 770 G B$ and the chaotic function x(t). 

Performing a similar procedure as presented above, by means of the solution ((£ — 6*1,771) of system (|2.10l) . 
one can obtain that the inequality 

Ikwc^o) - at <kM\ * 4r ( eL/(no+2)T - : ) 

holds for all t £ [01,6*1 + (710 + 2)T]. On the other hand, the existence of a number £2(771) G [0, 2T] such that 
((£2(771), ?7i) G 5 is provable. 
Therefore we have 

\\r) x (t)(2n T + ti(r) ) + t 2 (r)i),r)o) ~ C(n T + £2(771), m)\\ < e , 

and hence 772 = r)x(t)(^2,r)o) € -B5, where 6» 2 = 2n Q T + £1(770) + £2(771). 

One can continue in the same manner to constitute a sequence {tj} , which depends both on the initial 
condition 770 £ B$ and the chaotic function x(t), and satisfies < t j < 2T, j > 1, such that for each integer 
7 > we have 77^ = r] x ^(9i, 770) S -Ba and 

|k (t )(£, - C(t - 0i, < ^ ( e M"o + 2)T - l) , (3.29) 

where the sequence {0i} is defined through the equation 

i 

6, = in T + ^2t j: i>l, (3.30) 
j'=i 

and 6*o = 0. We emphasize that for any i > 1, it is true that 6>; S [inoT, i(no + 2)T] and 6*i — flj-i = n$T + ti < 
(no + 2)T. The procedure of the proof for £ S [6*i, 6^+1] is illustrated in Figure [3] 

In the remaining part of the proof, we will indicate boundedness of the function r] x (t) (£, 770) — p{t), which 
implies boundedness of r) x ^(t,r/o). We note that 77 x (i)(£, 770) — p(t) is a solution of system (|2.15|) . satisfying 
the initial condition r) x /£\ (0, 770) — p(0) = rjg. 

Using the relations 

at-o j ,v j ) = m+ [ nc(s-e jt ru))ds 

and 

p(t)=p(9 j )+ [ f( P ( S ))ds 

on the time interval [0j , 9j + 2T] , together with Gronwall-Bcllman Lemma, we achieve the inequality 

\\at-9 J , Vj )-p(t)\\<(6 + p)e^ T . 
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Figure 3: The representational visualization of the proof of Theorem 13. II The trajectory in red color shows 
the solution r] x (t) (t,r]o) of system p. lip , while the trajectory in blue color represents the solution £(t — 6i,rn) 
of system (|2.10p . where the sequence is defined through formula (|3.30[) and m = 77 x (t)(0jj 770). The 
presented illustration covers the way of the 2— dimensional case of the proof on the time interval [9i,6i+i], 
for an arbitrary nonnegativc integer i. At the moment t = the solution £(i — Oi,rji) belongs to the ball 
Bis, and on the other hand, rj i+1 is inside the ball B$. 

Hereby, one can see that 

HCfe+ife)>%)ll < Kitj+xivjUj) -p(t j+ i(vj))\\ + lb(W^))ll 

< P +{5 + p)e 2L i T . 

Since the point C(tj+i(r)j),rij) is on the surface S, according to inequality (|3.26[) . for all t > 0, we have 

\\C(t + t J+ i(vj),Vj)-p(t)\\ = ||C(t,Cfe+ife),%))-K*)ll 

<4A- 1 ||p- 1 (0)||||C(t, + i(^), ?7j )||e- Q */ 2 
< 4A X ||P _1 (0)|| [p+(S + P ) e 2L f T ] e~ at ' 2 
and therefore one attains for all i > that 

IICC* + tj+i(Vj),Vj) - P(t)\\ < iK i ll P_1 (°)|| [P + ( S + P) e2LfT ] ■ 
Thus, for t e [9j, 9j+i], we have 

\\iix{At-,m) - p{t)\\ < \\vx(t)(t,vo) - C{t - 9j,Vj)\\ + IICC* - Qj,Vj) -p(t)\\ 

where H(8,p) is the maximum of the numbers (8 + p)e 2L f T and AK\ ||P _1 (0)|| [p + (8 + p)e 2LfT ] . It is 
worth saying that H(S,p) — > as 8 — > and p — > 0, and ||?y s (t)(t, ?7o) — p(t)\\ can be made arbitrarily small 
by suitable choices of p, 8 and p. 

Consequently, any solution r] x ^(t, 770) of system (|2.1ip satisfying the condition 770 G Bg is bounded on 
R+. 

Proof of the theorem is finalized. □ 

Now, we fix values of p and 8 as specified in Theorem 13. 11 and use these values in the next section, where 
we will discuss the extension problem of sensitivity. 

4 Entrainment of sensitivity 

In the present section, extension of sensitivity feature through system (12. lip will be handled. Wc shall begin 
with the meaning of the aforementioned property for systems (|2.9p and (12. lip , and then state the main 
theorem of the present section. 
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We say that system (|2.9[) is sensitive if there exist positive real numbers eo and A such that for an arbitrary 
positive real number So and for each chaotic solution x(t) of system (|2.9p . there exist a chaotic solution x(t) 
of the same system and an interval J C R+, with length not less than A, such that ||a;(0) — x(0)|| < 5o and 
\\x(t) - x(t)\\ > e , for all t e J. 

Now, suppose that S is a positive number as indicated in Theorem l3.ll In a similar way, the system (|2.1ip 
is called sensitive if there exist positive real numbers e\ and A such that for an arbitrary positive number 
5\ < S and for each solution T] x (t)(t, rjo), where x(t) is a chaotic function and rjo € B$, there exist an interval 
J 1 C J, with length not less than A, and a solution 7fe(t)(t, 771), where x(t) is a chaotic function and rji £ B$, 

such that 1 1 770 - m\\ < s i aild \\Vx(t)(t,Vo) ~ > £i, for all t e J 1 . 

Through the next theorem, the extension of sensitivity feature from system (|2.9[) to system (|2.11j) will be 
mentioned. 

In addition, we shall say that system (|2.9[) + (|2.11[) is sensitive if systems (|2.9p and (|2.11[) are both 
sensitive. This description is a natural one since, otherwise, the inequality \\x(t) — x(t)\\ > eo implies that 
|| (x(t), r] x (t)(t)) ~ Vx(t) (t)) || > e o in the same interval of time, which already signifies unpredictability 

of system (|2.9|) + (|2.11j) . In our theory of cntrainment of chaos the crucial idea is not only the extension of 
sensitivity through system (|2.9p + (|2.11[) . but also through system (|2.11[) individually. For this reason, one 
should understand sensitivity of system (|2 . 9[) + (|2.11[) as a property which is equivalent to the sensitivity of 
system (|2.11|) . 

Theorem 4.1 If system (|2.9p is sensitive then the same is true for system (|2.1ip . 

Proof. Fix an arbitrary positive number 5i < 5 and let rj x / t \(t,r)o) be a given solution of system (|2.1ip . 
where x(t) is a chaotic function and 770 G B$. Since system (|2.9p is sensitive, there exist positive real numbers 
eo and A such that for any fixed number Sq > the inequalities \\x(0) — x(0)\\ < 5o and \\x(t) — x(t)\\ > eo, 
t e J, hold for some chaotic solution x(t) of system (|2.9p and for some interval J C M+ , whose length is not 
less than A. 

Now, let us fix 771 € B$ such that 1 1 770 — r\\ | < S\. In the proof our aim is to determine positive real numbers 
ei and A such that for some interval J 1 C J with length A, the inequality H^^ft) (i 5 5?o) — >fe(t)(^> , 7i) > e ii 
holds for alH G J 1 . 

/ 9r{x) \ 



Suppose that ^(x) 



where each ^ , 1 < j < n, is a real valued function. 



\ 9n(x) J 

Since for each chaotic solution x(t) of system p.9p the function x'(t) is inside the tube with radius Mp, one 
can conclude that the collection of chaotic solutions of system (|2.9p constitute an equicontinuous family on R + . 
Making use of the uniform continuity of the function ~g : M. m x M. m — > M. n , defined as g(xi, 22) = g{xi) — 3(^2), 
on the compact region 

S = {(x u x 2 ) e M m x K m I < Af, ||x a || < M} , 

together with the equicontinuity of the collection of chaotic solutions of system (|2.9p . one can verify that 
the set consisting of the elements of the form gj(x(t)) — gj(x(t)), 1 < j < n, where x{t) and x(t) are chaotic 
solutions of system (|2.9p , is an equicontinuous family on M. + . 

Therefore, there exists a positive real number r < A, independent of the functions x(t) and x(t), such 
that for any £i,f2 € R+ with |tj — ^| < t, the inequality 

I (.9, (x^)) - 9j (x( tl ))) - ( 9j (x(t 2 )) - 9j (x(t 2 )))\ < (4.32) 

holds for all 1 < j < n. 

Condition (A2) implies for alH € J that ||g(x(i)) — g(x(i))|| > L g \\x(t) — x(t)\\ . Therefore, for each t € J, 
there exists an integer jo, 1 < jo < n, which possibly depends on t, such that 

\g JO (x(t))~g JO (x(t))\>^Mt)-m\\- 
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Otherwise, if there exists s <G J such that for all 1 < j < n, the inequality 

\9j{x(s))-gj(x(s))\ < -^-\\x(s)-x{s)\\ 
holds, then one encounters with a contradiction since 

n 

\\g(x(s)) - g(x(s))\\ < \9Ms)) - 9j(x(s))\ < L g \\x(s) - x(s)\\ . 

3 = 1 

Now, let So be the midpoint of the interval J and 9 = sq — ^. One can find an integer jo = jo{so), 
1 < jo < n, such that 



\9h (x(s )) - g jo (x{s ))\ > \\x(s ) - x(s )\\ > 



L g e 



(4.33) 



On the other hand, making use of the inequality (|4.32[) , for all t £ [9, 9 + r] we have 

\9j 0(s )) ~ 9j (x(s ))\ - \g JO (x(t)) - g JO (x(t))\ 

< \ (9j ~ 9 j (x{t))) - (g jo (x(s )) - g jo (x(s )))\ 

„ L a e o 



2n 



and therefore, by means of (|4.33[) . we achieve that the inequality 

\9j ( x (t)) ~ 9j (»(*))! > Iflio 0(so)) ~ gjo (x~(s ))\ 



L a e o > L a ( o 
2n 2n 



(4.34) 



is valid for all f £ [6, 9 + r] . 

Making use of the mean value theorem for integrals, one can find numbers s\, s 2 , . . . , s n € [9,9 + r] such 



that 



[g(x(s)) -g(x(s))]ds 



[gi(x(s)) -5i(x(s))]ds 



1 + T 



[g2(x(s)) -g 2 (x(s))]ds 



0+T 



[g n {x(s)) - g n {x(s))]ds j 

I r[gi(x(si)) -gi(x(si))] ^ 
r [fl> 2 (x(s 2 )) - t/ 2 (x(s 2 ))] 



\ r [3„(x(s„)) - g n {x(s n ))] j 
Thus, using the inequality (|4.34|) . we attain that 

9 + T 



[g(x{s)) -g(x(s))]ds 



> t \g jo (x(s jo )) - 9j (x(s jo ))\ > 



2n 



(4.35) 



For t € [9, 6 + t], the solutions f? x (i)(^) 7 ?o) an d 7 te(t)(i J ??i) satisfy the integral equations 



Vx(t){t,Vo) = Vx(t){9,Vo) + / f(Vx(t){s,Vo))ds + fj,g(x(s))ds 
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and 

Vx(t)(t,m) =%(t)(Mi) + / /('fe(t)( s ,'?i)) ds + / vg(x(s))ds, 

Je Je 

respectively, and herewith the equation 

Vx(t)(t,Vo) ~ fte(t)(*,f7i) = {Vx(t){0,Vo) - Jte(t)(0,»7i)) 
+ / [/(%(i)(s,??o)) -/(%(t)(s,?7i))] «f« 

+ / - g{x{s))]ds 

Je 

is achieved. Hence, we have the inequality 

i-e+r 

||»7»(t)(^ + T 5 »7o) -Vx(t)(0 + T)'7i)|| > M / [#(>(«))- ff(z(s))]ds 

Je 

- ||»7x(t)(^»7o) -»te(t)(^m)|| - / L / ||%(t)(s,??o) - %(*)(«, > 1 1 ds. 

Je 

Now, assume that 

ImI TL g e 



(4.36) 



max ||77 K ( t )(t,?7o) - %(*)(*, t?i)| 



< 



In this case, one arrives at a contradiction since, by means of the inequalities (|4.35[) and (|4.36[) . we have 

1 6™ e+ T ] 1 1 ^ (*) ) ~ ^(t) ^! ) 1 1 - 1 1 ^ (t) (0 + r > ? ?0 ) - %(t) + T, J7l ) 1 1 

> ^^ Z/g£ ° - (1+rX/) max l|r? x( t)(i, r? ) - »?i)|| 

2n t£[e,e+r] " w w " 

> 



2n v -^2n(2 + TL/) 



2n V 2 + tL I 

(Ml ^igCO 



2n(2 + ri / ) 

Therefore, we have 

11 /- \ /nii lMl r ioeo 

t m^ll^^r/o)-^)^^)! > 2 „ (2 + t£/) - 

Suppose that the function ||?7 x (t)(t, %) — ^te(t)(^ 7 ?i)|| takes its maximum on the interval [6,0 + t] at the 
point £, that is, 

max 1 1 7?o) - »te(f)(*,57i)|| = \\Vx(t)(£,,Vo) ~ 'te(t)(£,*7i)|| > 

te[8,8+r] 1 v " 

for some # < £ < # + t. 

For i G [0, + t], by favour of the integral equations 

Vx(t)(t,Vo) = Vx(t)(£,Vo) + / f(Vx(t)(s,Vo))ds+ fig{x(s))ds, 
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and 



we obtain 



Define 



and let 



+ J [f{Vx(t)(s,m)) ~ f(Vx-(t){s,Vi))] ds + ^9(x(s)) - 5(S(s))]ds. 



2' 8n(K L f + M g | j u|)(2 + rL/) 



£-r\ if£># + f 



We note that the interval J 1 = 9 1 + r 1 ] is a subset of [6, 9 + t] and hence a subset of J. 
For t E J 1 , we have 

||»?s(i)(<,»7o)-»te(i)(*,»7i)|| > ||^(t)(^%) ~ %(t)(C,»7i)|| 



> 



> 



Lf ||^(t)(s,77o) - %( t )(s,?7i)|| d 

-2A" i/r 1 -2A/ 9 |/i|r 1 
-2T 1 {K Lf + M a \n\) 



-\iA 


I 







\g(x(s)) -g(x(s))\\ ds 



2n(2 + ri/) 

ImI T^ 3 eo 
2n(2 + ri/) 

H i~L g e 



4n(2 + rL/)' 
Consequently, we achieve for t E J 1 that 

||%(t)(*,??o) - nx(.t)(t,Vi)\\ > ei, (4.37) 

and the length t 1 of the interval J 1 does not depend on the chaotic functions x(t) 



where e\ 



\lA rL g e 



4n(2 + TLf) 
and x(t). 

The theorem is proved. □ 

Thus the important property of sensitivity is proved. This property can be considered as the unique 
ingredient of chaos for a set of bounded solutions [22j [32] . 

The formula (|3.29|) provides us a support that chaotic solutions are near solutions of non-perturbed 
cyclically behaved motions and consequently, for sufficiently small , the system admits chaotic cycles. In 
the same time, in general, we do not request necessarily appearance of cyclical properties for chaotic solutions. 
That is, if is not sufficiently small, then chaotic solutions may not be cyclical. So, our results provide this 
possibility also. 

Now, our aim is to exemplify our theoretical results through the illustration of 2— dimensional projections 
of the chaotic attractor and the Poincare section of a sample 4— dimensional system which is in the form 
of system (|2.9|) + (|2.11j) . We will also verify through simulations the existence of the sensitivity feature in 
accordance with Theorem 14.11 

For our purposes, we shall take into account the following 2— dimensional system 



u[ = aui — it 2 — Ui(u\ + itj) 
u' 2 = U\ + au 2 — 1*2(1*1 + 1*2), 



(4.38) 
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which is in the form of system (|2.10p , where a is a positive real number and 



One can verify in this case that p(t) = ( ^/^g^^ J is a periodic solution of system (|4.38[) . Evaluating 



A(t) = f'(p(t)), we achieve that 



/ -2acos 2 i -l-asin(2i) \ 
A(t)=[ 2 . (4.40) 

\ 1 -a-sin(2f) -2a sin 2 1 J 

Since the variational system (|2.13j) . where the matrix function A(t) is given by formula (|4.40p . admits 
the number pi = 1 as a multiplier, Lemma 7.3 |33] implies that the second multiplier is given by p 2 = 

cxp ( J Q 27r trA(s)ds S j = e~ 4lrQ . Thus, the second multiplier of the variational system corresponding to system 

(|4.38|) is in modulus less than one, and consequently the periodic solution p(t) of this system is asymptotically 
orbitally stable according to Andronov-Witt Theorem. 

We take into consideration the chaotic DufHng's oscillator [25], as the generator of chaos, presented by 
the differential equation 

x" + 0.05a:' + x 3 = 7.5 cos t. (4.41) 
Defining the variables x\ = x and x 2 = x' , equation (|4.41[) can be rewritten as a system in the following form 

x' 2 = -0.05a; 2 - x\ + 7.5 cost. v ' 

In the remaining part of the example, we make use of the value a = 9 in system (|4.38j) . and set up a 
unidirectional coupling between the systems (|4.42p and (|4.38|) to achieve the following 4— dimensional system 



x\ = x 2 , 

x' 2 = —0.05x2 — xf + 7.5 cost, , . 

x' 3 = 9^3 — X4 — Xy,{x\ + xl) + 0.5X1, 
x 'i = x 3 + 9^4 — £4(23 + xl) + 3.6X2. 



We suppose that system (|4.43|) admits a chaotic attractor in the 4— dimensional phase space. In Figure 
21 we visualize the 2— dimensional projections on the x\ — x 2 and X3 — X4 planes of the trajectory of system 
with initial data xi(0) = 3.05, x 2 (0) = 4.153, x 3 (0) = 2.8, 2:4(0) = 0.5. The picture presented in Figure 
IH (a) is in fact the chaotic attractor of the system (|4.42l) and Figure HI (6) represents a chaotic attractor in a 
neighborhood of the limit cycle of system (|4.38j) . with a = 9. It is observable that the attractor around the 
limit cycle shown in Figure HI (b) repeated the chaotic structure of the Ueda's attractor presented in Figure 
II (o). 

Figure [5] pictures the 2— dimensional projections of the whole Poincare section inside the 4— dimensional 
phase space of system (|4.43p . which is obtained by marking the trajectory corresponding to initial data 
xi(0) = 2,2:2(0) = 3,2:3(0) = 3,2:4(0) = stroboscopically at times that are integer multiples of 27T. Figure 
(a) represents the projection of the Poincare section on the 2:4 — 2:2 plane, and we note that this projection 
is, in fact, the strange attractor of system (|4.42p . On the other hand, the projection on the 2:3 — 2:4, plane 
presented in Figure (b) is a strange attractor around the limit cycle of system (|4.38p . It is apparent that 
the attractor indicated in Figure El (b) repeated the structure of the attractor shown in Figure EJ (a) . The 
simulation results presented in Figure [4] and Figure [5] indicate the concept of entrainmcnt of chaos. 

Now, to show through simulations the extension of sensitivity, we consider two initially close solutions 
of system ([£35]) . one with the initial data xi(0) = 3.07, x 2 (0) = 4.18, x 3 (0) = 1.57, x 4 (0) = -0.25, which is 
presented in blue color, and another with the initial data xi(0) = 3.22, 2:2(0) = 4.14, x 3 (0) = 1.35, X4(0) = 
—0.22, which is pictured in red. In Figure [SI we present 2— dimensional projections of these trajectories on 
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a b 




Figure 4: The 2— dimensional projections of the chaotic attractor of system (|4.43[) . (a) Projection on the 
x\ — X2 plane, (b) Projection on the X3 — x& plane. The theoretical results indicate that system (|4.43|) admits 
a chaotic attractor in the 4— dimensional phase space. The picture in (a) represents not only the projection 
of the attractor of system (|4.43p on the x\ — X2 plane, but also the chaotic attractor of system (|4.42|) . The 
chaotic attractor pictured in (b), which appears around the limit cycle of system (|4.38p with a = 9, is a 
manifestation of the entrainmcnt mechanism. 

the X\ — x-i and X3 — X4 planes. The picture in Figure [5J (a) shows the sensitivity feature of system (|4.42p . 
while picture in Figure [SI (b) represents the extension of this feature. 

We note that formula (|4.37[) implies that the strength of sensitivity of system (|2.1ip is proportional to the 
strength of the chaotic perturbation, ng(x), used in this system. Therefore, if one considers system (|2.1ip 
with weak perturbations, in that case, although the extension of sensitivity feature is guaranteed by Theorem 
14.11 it may not be visible in simulation results. But, according to formula (|3.29p . strong perturbations may 
diminish the cyclic behavior of the chaotic solutions. In this sense, the picture presented in Figure [HI (b) 
exhibits the extension of sensitivity, but does not indicate a cyclic behavior for the illustrated solution, 
according to the strength of the perturbation used in system (|4.43|) . 

5 Entrainment of period-doubling cascade 

Period-doubling route to chaos is nowadays a popular way of indication of chaos |34|-|37| which has applica- 
tions in mechanics, electronics, biology, and is most convenient for simulations |38|-|40j. 

We start this section by describing the chaos through period-doubling cascade for system (|2.9p and 
continue with its extension through system p. lip . 

To discuss the existence of chaos through period-doubling cascade, we consider the system 

x' = G(t,x,X), (5.44) 

where A is a parameter and the function G : M. + x R m x R — > R m satisfies the property that F(t,x) = 
G(t, x, Aoo), for all t e R + , x E R m and for some finite value Aoo of the parameter A, which will be explained 
below. 

The following condition is needed throughout the section. 

(A4) There exists a positive real number u> such that the periodicity property G(t + lo,x, A) = G(t,x, A) 
holds, for all tel + ,i6 R m and A <G R. 

System (|2.9p is said to be chaotic through period-doubling cascade if there exist a natural number fc 
and a sequence of period-doubling bifurcation values , \j — > Xoo as j —> 00, such that for each natural 
number j as the parameter A increases or decreases through Xj, system (|5.44|) undergoes a period-doubling 
bifurcation and a periodic solution with period k^uj appears. As a consequence, at the parameter value 
A = Aqo, there exist infinitely many unstable periodic solutions of system (|5.44p . and hence of system (|2.9I) . 
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Figure 5: The 2— dimensional projections of the whole Poincare section inside the 4— dimensional phase space 
of system (|4.43[) . (a) Projection on the x\ — X2 plane, (b) Projection on the x$ — X4 plane. The image shown 
in (a) is, in fact, the Ueda's strange attractor obtained through the Poincare section of system (|4.42[) . The 
picture in (b) represents a strange attractor around the limit cycle of system (|4.38|) . Possibly one can call the 
strange attractor presented in (b) as the Ueda 's strange attractor near limit cycle. 

all lying inside the tube with radius M. In this case, system (|2.9|) admits periodic solutions of periods 
k uj, 2k u, 4:k oj, 8k oj, ■ ■ ■ . For details one can see [Ml |4"TI|4"2"] . 

In a similar way, we say that system (|2.11[) admits the chaos through period-doubling cascade if system 
(|2.9[) is chaotic through the same route and for each k^uj— periodic solution Xj(t) of system (|2.9[) . there 
exists a periodic solution (t) (t, f?o) of the system 



with the same period such that all such periodic solutions lie in a bounded region. Accordingly, system (|2.11|) 
possesses infinitely many periodic solutions with periods ko2 1 uj for each natural number j such that all of 
these periodic solutions are unstable according to Theorem 14. II 

In the case that system (|2.11| is chaotic through period-doubling cascade as described above, we say that 
system (|2.9|) + (|2. 1 1|) admits the chaos through the same route and extension of this type of chaos is provided. 

If system (|2.11|) is chaotic through period-doubling cascade, in the same way as system (|5.44j) . then 
system (|2.11[) also undergoes period-doubling bifurcations as the parameter A increases or decreases through 
the values Xj. That is, the sequence of bifurcation parameters is exactly the same for both systems. 
In the case system (|5.44[) , and consequently system (|2.9[) , obey the Feigenbaum's universal behavior (26] 
136] SU |43j |44], one can conclude that the same result holds, also, for system (|2.11[) . In other words, when 



achieved and this universal number is the same for systems (|2.9[) and (|2.11[) . and consequently for system 

mu) + dun). 

Now, suppose that x(t) is a po— periodic solution of system (|2.9p for some positive real number pq. 
Theorem (|3.ip indicates that a bounded solution of the system y' = f(y) + g(x(t)) exists, provided that 
is sufficiently small. In the case that the dimension of system (|2.1ip is n = 2, Massera's Theorem [15] 25] 
implies the existence of a po —periodic solution of the same system. On the other hand, by condition (^42) 
the converse is also true. That is, if there exists a po— periodic solution of the system y' = f(y) + g(x(t)), 
then x(t) is necessarily po~ periodic. Therefore, if system (|2.9|) admits infinitely many periodic solutions 
with periods k&uj for each natural number j, then the same is true for system (|2. 1 1|) provided that n = 2. 
Additionally, no periodic solutions with any other period exist for this system. 

Next, we continue with an example which provides extension of chaos through period-doubling cascade. 
In paper [57], it is mentioned that the Duffing's equation 



v = f(y) + pg(xj(t)) 



(5.45) 



lim,- 




is evaluated, the universal constant known as the Feigenbaum number 4.6692016 ... is 



x" + 0.3.x' + x 3 = Acosi, 



(5.46) 
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a b 




X 1 X 3 

Figure 6: Extension of sensitivity in system (|4.43p . (a) Projection on the x\ — X2 plane, (b) Projection on 
the X3 — X4 plane. The picture in (o), in fact, represents two initially nearby trajectories corresponding to 
system (|4.42[) , which eventually diverge. The extension of sensitivity is observable in (b) such that the blue 
and red trajectories are initially close to each other and then separated. It is noteworthy that since the 
chaotic perturbation in system (|4.43|) is not weak enough as indicated in Theorem I3.ll it is visible in picture 
(b) that the cyclic behavior is not achieved. The reason for the usage of a strong perturbation is to increase 
the strength of sensitivity, since otherwise, it is not able to visualize extension of this feature. 



where A is a parameter, displays period-doubling bifurcations and admits the chaos through period-doubling 
cascade at the parameter value A = Aoo = 40. Making use of the new variables x\ — x and x 2 = x', we can 
rewrite equation (|5.46|) as a system in the form 



X2, 



-0.3^2 — xf + A cost. 



(5.47) 



To illustrate entrainment of chaos through period-doubling cascade, we shall make use of the parameter 
value A = A^ in (|5.47[) and combine with system (|4.38|l to constitute the 4— dimensional system 



x 2 , 



-0.3x 2 -x\+ 40 cost, 

axz — X4 — x 3 (x§ + X4) + fixi, 

x 3 + ax4 — X4 {x\ + x\) + [1X2, 



where /i is a non-zero real number and a is a positive number. 

We shall make use of the technique of Lyapunov functions to identify an annular region G C 
contains infinitely many periodic solutions of system 



y'i 
y'2 



ayi - 2/2 
Vi + ay 2 



yi \y\ 



yi \v\ 



2 
2 

yi 



fJLXi, 



(5.48) 



which 



(5.49) 



which is in the form of system (|2.11[) . and this discussion will imply the existence of chaos through period- 
doubling cascade in that region. In such a case, the 2— dimensional projection on the X3 — X4 plane of the 

chaotic attractor of system (|5.48[) will appear around the orbit 7 of the periodic solution p(t) = 

of system fOS]). 

Consider the Lyapunov function V(x3,X4) 



/asint 



c? + x\. In this case one has 



%^gl (^3, X4) = 2x3X3 + 2x4X4 
= 2x3 \0tX3 — X4 




19 



The inequality x\ + x\ > \fa + ri, where t\ is a positive number, implies that 



x\ + x\ \x\ + x\ - a) > rf + 3rf^/a + 1t\Ol. 



Similarly, if < r 2 < sJx\-\- x\ < s/a — r$ , where r 3 < -^/a, then we have 



' x\ + x\ {x\ + x\ - a) < r 2 r\ - 2r 2 r 3 y/a. 
Since the chaotic attractor of system (|5.47p satisfies \xi \ < 6 and |cc2 1 < 15, we attain that 

xix 3 x 2 X4 



\J x\ + x\ \J x\ + x\ ^ 
Therefore, if is sufficiently small such that 



<21|/x|. 



ImI < 21 max { r l^ r Vi 



(5.50) 



then one can verify that ^ j fe^gt (^3; X4) < for \J x\ + x\ > \fa + r\ and V fe ^ (0:3, 0:4) > for < r 2 < 
\J x\ J r x\ < y/a — 7'3 . Consequently, there are infinitely many periodic solutions inside the region 



G= I (0:3,0:4) El 2 : \foi — ?* 



3 < \l x\ + x\ < y/a + ri 



It is worth saying that the region G can be made arbitrarily narrow by small choices of the numbers n, r 3 
and |//| according to inequality (|5.50[) . 

Now, let us consider system (|5.48p with a = 0.002 and /j, = 0.008 and take into account a trajectory 
of this system corresponding to initial data xi(0) = 3.5, 0:2(0) = — 2,0:3(0) = 0.02 and 0:4(0) = 0.038. The 
2— dimensional projections of this trajectory on the x\ — x 2 and x 3 — X4 planes are presented in Figure [7J 
The picture in Figure [71 (a), illustrates the chaotic attractor of system (|5.47p and the picture in Figure [JJ (b), 
represents the chaos around limit cycle. 





Figure 7: 2— dimensional projections of the chaotic attractor of system (|5.48p . where a = 0.002 and pt = 0.008. 
(a) Projection on the x\ — x 2 plane, (b) Projection on the 0:3 — 0:4 plane. The picture in (a), in fact, represents 
the chaotic attractor of system (|5.47|) . The extension of chaos through period-doubling cascade is observable 
in picture (b) such that the projection on the 0:3 — X4 plane represents chaotic behavior. Appearance of the 
attractor in picture (b) around the limit cycle is a consequence of the entrainment process. 



Note that system (|5.48|) exhibits a symmetry under the transformation 

: (x±, x 2 , x 3 , x 4 , t) — !• {—x\, ~x 2 , —0:3, —0:4, t + vr). 
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Therefore one can conclude that the chaotic attractor in the 4— dimensional phase space of system (|5.48p is 
symmetric around the origin, and that is the reason for the symmetry of the projections on the x\ — X2 and 
X3 — X4 planes presented in Figure [7] 

In the next theorem, we will generalize our discussions about extension of chaos through period-doubling 
cascade in the case that system (|2.1ip is n— dimensional for arbitrary natural number n. Before we state the 
theorem, we introduce the following definition. 

We say that the solutions of the non-autonomous system 

v' = h(t,v), (5.51) 

where h : R x R 9 — > R q is a continuous function in all its arguments, are ultimately bounded for bound B 
if there exists a number B > such that for every solution v(t),v(to) = vq, of system (|5 . 5 1[) . there exists a 
number R > such that the inequality ||t>(i)|| < B holds for all t > to + R, where B is independent of the 
particular solution while R may depend on each solution. 

The proof of the next theorem can be verified using Theorem 15.8 [46] . 

Theorem 5.1 If system (|2.Q[) admits the chaos through period- doubling cascade and if there exists a positive 
number B such that for each chaotic solution x(t) of system (|2.9|) solutions of system y' = f(y) + jJ,g(x(t)) 
are ultimately bounded for bound B , then system (|2.1ip is chaotic in the same way. 

We again emphasize that in the case of extension of chaos through period-doubling cascade, the unstability 
of the infinite number of periodic solutions of system (|2.1ip is assured by Theorem 14. II On the other hand, 
one can see that conditions of Theorem (|5.ip are fulfilled according to our discussions through Lyapunov 
functions for system (|5.48p , and this also provides existence of chaos through period-doubling cascade in the 
system. 

6 Discussion 

We shall devote this section to discuss through simulations the problems of entrainment of chaos by toroidal 
attractors, entrainment in Chua's oscillators and chaos control problem. We start with the demonstration of 
chaos generation around tori. 

6.1 Entrainment of chaos by toroidal attractors 

In previous parts of the paper, we have discussed entrainment of chaos with limit cycles. Now, the question 
is whether a similar approach is possible around tori. In this part, we will consider, numerically, the problem 
of seizure of chaos by tori. 

Let us consider the following 3— dimensional system [35J |H] 

ui = (A - 3)ui - 0.25it 2 + ui (u 3 + 0.2(1 - uj)) , 

u' 2 = 0.25ui + (A - 3)u 2 + u 2 (m + 0.2(1 - u\)) , (6.52) 
u' 3 = Xu 3 - (ul + ul + ul), 

where A is a parameter. 

For small and positive values of the parameter A, system (|6.52|) admits an asymptotically stable equilibrium 
point, with a positive u 3 coordinate near the origin, and at A ~ 1.68, the equilibrium point loses its stability 
and an hyperbolic, asymptotically orbitally stable limit cycle takes place. At the parameter value A = 2, 
the periodic orbit is still asymptotically orbitally stable, but not hyperbolic. For A > 2, the limit cycle is 
no longer stable and an attracting invariant torus takes place near the periodic orbit. With the increasing 
values of A, the invariant torus grows rapidly [48] . 

Figure E] illustrates the trajectory of system (|6.52p with the parameter value A = 2.003, corresponding to 
the initial data ui(0) = —0.0983, 1*2(0) = 0.9004 and 143(0) = 0.6908. It is seen in the figure that the motion 
of the trajectory is around a torus. 
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Figure 8: The attracting invariant torus of system (|6.52[) with A = 2.003. 



To achieve chaotic motions around torus, we perturb system (|6.52[) through the solutions of the chaotic 
Lorenz system [50] 

x[ = — 10xi + 10x2, 

x' 2 = —X1X3 + 28:ri — x 2 , (6.53) 
x' 3 = xix 2 - (8/3)2:3, 

and constitute the following 6— dimensional system 

x[ = — 10xi + 10x2, 

x' 2 = —X1X3 + 28xi — x 2 , 

x' 3 = X1X2 - (8/3)x 3 , 

x' A = (A - 3)x 4 - 0.25x5 + x 4 (x 6 + 0.2(1 - x 2 6 )) ( . 

4 = 0.25.T4 + (A - 3)x 5 + x 5 (x 6 + 0.2(1 - x 2 6 )) 
+0.004x2, 

x' 6 = Ax 6 - (x| + x\ + x%) + 0.002x3, 

where A — 2.003 once again. Let us consider the solution of system (|6.54|) corresponding to the initial data 
xi(0) = -6.7453, x 2 (0) = 0.3435, x 3 (0) = 32.7629, x 4 (0) = 0.0793, x 5 (0) = -1.1761 and x 6 (0) = 0.9449. 
Figure |^] shows the projection of the considered trajectory on the X4 — X5 — Xq space. It is observable in 
Figure IH1 that the motion is disposed to behave both chaotically and around a torus, simultaneously. 



6.2 Entrainment in Chua's oscillators 

We shall continue our investigation by presenting a simulation result for entrainment of chaos obtained in 
Chua's oscillators. The dimensionlcss form of a Chua's oscillator given by the system 

x[ = ka [x 2 - xi - ^(xi)], 

x' 2 = k(xi - x 2 + x 3 ), , . 

•1) a \ (6.55) 

x' 3 = k {-P0X2 -70X3), 

ij>(x) = bx + 0.5(a - b) (\x + 1| + \x - 1|) , 
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Figure 9: The seizure of chaotic motion by a toroidal attractor. In the figure, the chaotic behavior around 
the invariant torus of system (|6.52p with A = 2.003 is visualized. The picture is obtained by illustrating 
the 3— dimensional projection of the chaotic trajectory of system ()6.54|) . with initial data xi(0) = —6.7453, 
^(O) = 0.3435, 0:3(0) = 32.7629, x 4 (0) = 0.0793, x 5 (0) = -1.1761 and a: 6 (0) = 0.9449, on the x 4 - x 5 - x 6 
space. The appearance of chaotic motion around a torus represents the entrainmcnt process emerging in 
system (|6.54j) . 

where cto, Po, 7o ; a , b and k are constants. 

In paper [51], it is indicated that system (|6.55|) with the coefficients ao = 5 75 ;/3o = 7.8351, 70 = 
1.38166392 Q = _ 1 8 4 5 g ; 5 = _o.86604 and k = 1 admits a stable equilibrium. 

We revealed in Section [5] that the system (|5.48[) . with the coefficients a = 0.002 and fj, = 0.008, admits 
the chaos through period-doubling cascade and the 2— dimensional projection on the X3 — X4 plane of the 
chaotic attractor of this system appears near a limit cycle. Now, we attach two consecutive Chua's oscillators 
in dimensionless form, which admit asymptotically stable equilibrium points, to constitute the following 
10— dimensional system 

x[ = x 2 , 

x' 2 = — 0.3a; 2 — x\ + 40, cost, 
x' 3 = 0.002x3 - x 4 - x 3 (xl + xf) + 0.008a;i, 
x'i = x 3 + 0.002x4 - x 4 (X3 + xl) + O.OO8X2, 
x' 5 = fn^[x 6 - 0.13396x5 + 0.48993(|x 5 + 1| + \x 5 
Xg = x 5 — x 6 + x 7 + 2x 4 , 

x' 7 = -7.8351X 6 - 1 -38166392 X7 + ^ 

x' s = f^[x 9 - 0.13396x 8 + 0.48993(|x 8 + 1| + |x 8 

Xg = X 8 - Xg + X10 + 0.3X7, 

x' 10 = -7.8351x9 - 1 - 381 1 f 392 xi + 0.8x 6 . 

Since the subsystem with coordinates (xi,X2) has chaos through period-doubling cascade, the subsystem 
(x3,X4) admits the seized chaos by its limit cycle, which is guaranteed by our theoretical discussions. To 
illustrate an application of our results, we make use of system (X3, X4) as a source of chaos for the Chua circuit, 
presented by the subsystem (xs,X6,X7). In addition, the coordinates xs,X6 and X7 are used to perturb the 
next Chua's oscillator corresponding to the last three equations in system (|6.56[) . Following the results of 
paper [14], we have to observe chaotic behavior in both of the Chua's oscillators. 



1|)]+ 0.5x3, 
1|)]+ 0.8x 5 , 



(6.56) 
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We consider a trajectory of system (|6.56p with initial data xi(0) = 3.5, 2:2(0) = —2, x 3 (0) = 0.02, 
as 4 (0) = 0.038, 515(0) = -8.016, x 6 (0) = -0.084, x 7 {0) = 7.792, x 8 (0) = -22.764, ar 9 (0) = -0.281 and 
xio(0) = 20.167, and visualize its 3— dimensional projections on the X5 — xq — X7 and x$ — xg — X10 spaces 
in Figure [TO] We note that the projections on the x\ — X2 and X3 — X4 planes will give the same attractors 
presented in Figure [3 (a) and (6), respectively. The pictures presented in Figure [TOl (a) and (b) indicate that 
the chaotic Chua's attractors appear around limit cycles and this is a manifestation of entrainment of chaos. 
Furthermore, it is seen in Figure [TU] that the shapes of these attractors resemble the spiral Chua's attractor, 
which takes place in the case of a period-doubling cascade in Chua systems [5T1 152"] . 

a b 




Figure 10: 3— dimensional projections of the chaotic attractor of system (|6.56[) . (a) Projection on the X5 — 
xq — X7 space, (b) Projection on the x$ — xg — x\q space. The pictures presented in (a) and (b) both reveal 
that the Chua's oscillators possess chaotic behavior such that the chaotic attractors appear near limit cycles, 
and this is an indicator of entrainment of chaos. 



6.3 Controlling chaos 

The chaos obtained through period-doubling cascade admits infinitely many periodic solutions which are all 
unstable and Pyragas control method |53] is an effective instrument to stabilize them. One can say that it 
plays also an important role as it is very convenient to indicate through simulations the periodic solutions, 
which are invisible in the set of irregular motions. 

In this subsection, we will indicate a procedure by favour of an example for the stabilization of unstable 
periodic solutions of system (|2.9|) + (|2.1ip . The Pyragas control method will be mainly used in this procedure 
and therefore we start by explaining the method briefly. Pyragas, in 1992, proposed a delayed feedback 
control method for the stabilization of unstable periodic orbits of a chaotic system. In this method, one 
considers a system of the form 

x' = H(x,q), (6.57) 

where q = q(t) is an externally controllable parameter and for q = it is assumed that the system (|6.57|) is 
in the chaotic state of interest, whose periodic orbits are to be stabilized [36] [53] [54] [55] [56]. According to 
Pyragas method, an unstable periodic solution with period To of system (|6.57|) with q = 0, can be stabilized by 
the control law q(t) = C [s (t — tq) — s(t)] , where the parameter C represents the strength of the perturbation 
and s(t) = g [x(t)] is a scalar signal given by some function of the state of the system. 

It is indicated in [54] that in order to apply the Pyragas control method to the chaotic Duffing oscillator 
given by the system 

X l = X 2, / fi ro\ 

x' 2 = -0.10x2 + 0.5xi (1 - x?) + 0.24 sint, y ' 
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one can construct the corresponding control system 



1 — z 2, 

(6.59) 



z' 2 = -0.10Z2 + 0.52! (1 



+0.24sin(z 3 ) + C [z 2 {t - r ) - z 2 {t)} , 
4 = 1, 

where g(t) = C [22 (t — To) — z%(t)\ is the control law and the less unstable 2tt— periodic solution can be 
stabilized by choosing the parameter values C = 0.36 and To = 2ir. 

Making use of system (|6.58|) together with (|4.38p . where a = 7, we set up the following system 

x' x = x 2 , 

x' 2 = -0.10x2 + 0.5ari(l - x\) + 0.24 sint, . . 

%3 = 7x 3 — xa — x%(x\ + xf) + 5a?i, 

•4 = #3 + 7^4 — x 4 (a;§ + xf) + 4(x 2 + x 2 ). 

According to our theoretical discussions, system (|6 . 60[) admits a chaotic attractor in the 4— dimensional 
phase space, and its 2— dimensional projection on the X3 — X4 plane appears near the limit cycle of system 
(|Qg|l , with a = 7. 

Our present purpose is to demonstrate numerically how to control chaos of system (|6.60p . We propose 
that if a periodic solution of the 2— dimensional subsystem (2:1, £2), inside system (|6.60[) . is stabilized, then 
the chaos of system (|6.60p is controlled. In other words, it is enough to control the chaos of system (|6.58|) . 
which is used as the source of the exogenous perturbation in system (|6.60[) . 

To apply the Pyragas method for controlling the chaos of system (|6.60[) . we constitute the system 

A = Z2, 

z' 2 = -0.102 2 + 0.5zi(l - zf) + 0.24siii2 3 + C [z 2 {t - 2tt) - z 2 (t)} , 

4 = 1- (6-61) 

z\ = 7zi — Z5 - Z4,(zl + zf) + 5zi, 

4 = z 4 + 72 5 - 2 5 (2| + zf) + 4(22 + 4), 

which is the control system corresponding to system f|6 . 60|) . 

We consider the solution of system (|6.61[) with initial data 21 (0) = 0.1, 2 2 (0) = —0.8, 2 3 (0) = 0, 24(0) = 
2.64 and 25(0) =0.1. The system evolves freely taking C — until t = 70, and at that moment the control is 
switched on by taking C = 0.36. At the moment t = 210, we switch off the control mechanism and start to 
use the value of the parameter C = once again. Figure [11] pictures the graphs of the 22 and 25 coordinates 
of the solution which reveals the control of chaos of system (|6.60[) . It is also observable that after switching 
off the control mechanism, the stabilized 2ir~ periodic solution of system (|6.60[) loses its stability and chaos 
emerges again. Similar pictures can be obtained for the other coordinates of system (|6.61[) . which are not 
just pictured here. 



7 Conclusion 

The ideas of Huygens about synchronization are carried a step forward in such a way that the entrainment of 
chaos by limit cycles is discussed. Our theoretical results can be effectively used in arbitrarily high dimensional 
systems, which possess asymptotically orbitally stable limit cycles. Examples of such systems can be found in 
mechanics, electronics, economics, neural sciences, chemistry and population dynamics [31 HI [51 15T1 [551 1591 ISO] . 
Through the method presented, one can obtain motions which behave cyclically and chaotically in the same 
time. 

In the paper, theoretical problems with rigorous proofs for the existence of bounded solutions and ex- 
tension of sensitivity, which we consider as the unique ingredient of chaos, as well as extension of chaos 
obtained through period-doubling cascade arc considered. Additionally, the illustrated simulations support 
the theoretical results. Entrainment of chaos by toroidal attractors and entrainment in Chua's oscillators 
have been observed numerically. We indicate the existence of a chaotic attractor in a sample system by 
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Figure 11: Application of the Pyragas control method to system (|6.60[) by means of the corresponding control 
system (|6.61| . (a) The graph of coordinate, (b) The graph of z$ coordinate. The picture in (a) reveals that 
the control mechanism stabilizes the less unstable 2tt— periodic solution of system (|6.58j) . and correspondingly 
the picture in (b) supports that the chaos of system (|6.60|) is controlled. The control mechanism starts at the 
time t = 70 and ends at t = 210. On the intervals where the control is switched off, the irregular behavior of 
the coordinates are observable. 



means of Lyapunov functions. Moreover, the existence of unstable periodic solutions is discussed through the 
Pyragas method [53] and simulations. 

Cyclical behavior in chaotic attractors have been widely observed in the literature. We can refer for this 
famous Rossler attractor and Chua's spiral attractor, and even in the classical Lorenz attractor one can see 
two-center cyclical behavior. Our results for the achievement of cyclical behavior are different than those, 
since "cycling" chaos is usually obtained through period-doubling cascade, while our irregular cyclic behavior 
is a result of already existing chaotic motions, which are not in general "cyclic" , but applied as exogeneous 
perturbations. Nevertheless, one can guess that the mechanism proposed in our article may be underneaths 
of some chaotic attractors, which have already been discussed in the literature. 

The results of the paper can be mimicked for the case when cycles are attracting for the time decreasing 
to — oo. Moreover, they can be extended by considering tori as attractors. Another theoretically challenging 
problem is to consider hyperbolic cycles and also the critical cases. 

Our results are useful to generate multidimensional chaos, exceptionally if one requests that the phe- 
nomenon should be rigorously approved [fJTJ . If one considers Hopf bifurcation as a reason for the limit 
cycle generation, we can formally compare our results with the results of Ruelle and Takens [63] on the 
appearance of turbulence through three successive bifurcations. In our case, we have a chaos obtained after 
less than three bifurcations and additionally, we use chaotic perturbations. 
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